The Naimark dilated PT— symmetric brachistochrone 
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The quantum mechanical brachistochrone system with PT— symmetric Hamiltonian is Naimark 
dilated and reinterpreted as subsystem of a Hermitian system in a higher-dimensional Hilbert 
space. This opens a way to a direct experimental implementation of the recently hypothesized 
"PT— symmetric ultra-fast brachistochrone regime of [C. M. Bender et al, Phys. Rev. Lett. 98, 
040403 (2007)] in an entangled two-spin system. 

PACS numbers: 03.65.Xp, 03.65.Ca, 03.65.Ud, 03.67.Lx 



Introduction The quantum brachistochrone problem 
consists in finding a Hamiltonian H which evolves a given 
initial state jV'/) into a given final state \iPf) in a minimal 
time T. Considering this problem for quantum mechan- 
ics with T'T— symmetric Hamiltonians (PTQM) Bender, 
Brody, Jones and Meister (BBJM) found the surprising 
result 1] that the minimal evolution (passage) time r-pq- 
was less than the minimal time Th required for the evolu- 
tion induced by a Hermitian Hamiltonian [2, 3] . It could 
be made even arbitrary small r-pq- ^ in a strongly non- 
Hermitian regime [l|, |jl. If this effect of a 'faster than 
Hermitian' evolution [l| were experimentally realizable 
it would open a way to ultra-fast quantum computing 
processes [5|. A problem still unsolved in p| concerned 
the switching mechanism between the PT— symmetric 
brachistochrone system and a conventional (von Neu- 
mann) quantum system necessary for an experimental 
implementation of the suggested ultra-fast quantum pro- 
cess. 

As shown by Mostafazadeh [6(, an equivalence map- 
ping [71] between PTQM in the sector of unbroken 
T^T— symmetry and conventional quantum mechanics 
(CQM) leaves the passage time of a brachistochrone in- 
variant Tfi = T-pq-. This implies that a vanishing passage 
time Tpq- — > in the T^T— symmetric system is necessar- 
ily connected with a vanishing distance between initial 
and finite state in the equivalent Hermitian system — an 
effect geometrically analyzed in [8] . In case of the Hermi- 
tian equivalent of the BBJM brachistochrone, initial and 
final states will nearly coincide (coincidence problem) so 
that the brachistochrone effect in such an interpretation 
would loose any physical relevance. 

In this Letter, we propose a realization of the BBJM 
brachistochrone [1] which resolves the switching problem 
between PTQM and CQM regimes [l|, avoids the co- 
incidence problem [O, l8( and which can be considered 
as starting point for a direct experimental implemen- 
tation. The key idea consists in a reinterpretation of 
the BBJM brachistochrone as T'T— symmetric subsystem 
of a larger CQM system living in a higher-dimensional 
Hilbert space. For this purpose we use a Naimark di- 
latation (extension) technique (9| as it is widely used in 



quantum information theory [5|. We will demonstrate 
that the resulting large system will have the structure 
of an entangled two-spin (two-qubit) system so that an 
experimental realization of the BBJM brachistochrone ef- 
fect should be feasible, e.g., in a suitably designed system 
of entangled polarized photons [lO| . 

Technically, the construction of the large Hermitian 
system will be accomplished by a three-step procedure: 
(i) by building a suitable positive operator valued mea- 
sure (POVM) [5J, y, [llj over the nonorthogonal eigen- 
statcs of the T^T— symmetric Hamiltonian H and its ad- 
joint H\ (ii) by Naimark dilating (extending) [9|] this 
POVM into an orthogonal projector set in the higher di- 
mensional Hilbert space and (iii) by constructing from it 
a corresponding Hermitian Hamiltonian H = H^ and a 
unitary evolution operator U(t) 



-itH 



BBJM brachistochrone The BBJM brachistochrone 
1] that we are going to Hermitianly dilate (extend) de- 
scribes the evolution from an initial state \ipi) to a final 
state IV'-f) governed by a T^T— symmetric Hamiltonian H 
whose structure is chosen in such a way that the time r 
required for the evolution becomes minimal. As shown in 
i2] such a minimal-passage-time solution follows a mini- 
mal geodesic in projective Hilbert space and it is there- 
fore located in the two-dimensional subspace 7^2 = C^ 
spanned by \ijji) and lV'i=")- In this 7^2 the PT— symmetric 
Hamiltonian H can be chosen as [l|, |j] 



H = Eoh + s 



i sin(Q;) 1 

1 — isin(a) 



Eo,seR, (1) 



where V — <Tx denotes the parity operator, T is the an- 
tilinear operator of time reflection and complex conjuga- 
tion I4], Eq denotes an irrelevant offset energy and s a 
general scaling factor of the matrix. (As usual, cTx, (Jy 
and az are Pauli matrices.) The angle a S (— 7r/2, 7r/2) 
characterizes the non-Hermiticity of the Hamiltonian: 
H{a — 0) is Hermitian, whereas in the limit a -^ ±7r/2 
the Hamiltonian H becomes strongly non-Hermitian and 



similar to a Jordan block, i.e. its eigenvectors 



\E+{a)) 
\E-{a)) 



„Ja/2 



y/2cos{a) 
-\/2cos(a) 



(2) 



and eigenvalues E± = Eo±s cos(a) =: i?o±wo/2 coalesce 
for fixed \s\ < oo [J|. The Hamiltonian is restricted to 
purely real eigenvalues, i.e. the parameter sector of exact 
■PT-Symmetry ^. The operator U(t) = e""^ of the 
non-unitary evolution induced by H has the explicit form 



Uit) 



-iEot 



cos(a) 



cos{y — a) — isin(y) 
—i sin(y) cos{y + a) 



(3) 



with y := ujot/2 (we set h = 1). In the BBJM-brachis- 
tochrone setup Jj this U{t) is used to evolve an initial 
state \tpi) = (1,0)-^ into a final state \^pp) = fip{0,l)'^, 
jjLp :— —16^*^°"^. The time r required for this evolution 
follows from the condition w = a + 7r/2 as t = ° / s 

^ ' scos(a) 

and tends for 



e-7r/2, e^+0 



(4) 



and fixed scos(a) = ti;o/2 to zero: t ^ 0. In this way 
the evolution from \'tpj) to the orthogonal \4'f) induced 
by the T'T— symmetric Hamiltonian H with eigenstates 
of fixed energy difference -E+ — E- = ujq appears faster 
than an evolution between these states induced by any 
Hermitian Hamiltonian with the same energy difference 
loq between its eigenstates. This is due to the fact that 
the evolution time r between orthogonal states in Hermi- 
tian system has to be larger than the Anandan-Aharonov 
lower bound t > Th = tt/wq [3,|13|. 

Before we embed the BBJM-brachistochrone into a 
larger Hermitian model we briefly collect the required 
setup information. The eigenvectors ^ oiH are normal- 
ized with regard to the VT inner product (u, v) = VTu-v 
as {E±,E±) = ±1, {E±,E^) = 12] and for a ^ 
they are nonorthogonal with regard to the standard in- 
ner product in the Hilbert space 'H2 = C^: {E±\Ezp) ^ 0. 
We supplement them via H'^{a) — H{—a) with the eigen- 
vectors \E+{—a)), \E-{—a)) of the adjoint operator H'^ 
and arrange them as columns in the matrices 



* 



:|i?+(a)),|i?_(a))], ^■.= [\E+[-a))AE.(-a))] 



(5) 
With E := diag (£'+, E^) the eigenvalue problems for H 

and ijt take then the compact matrix form 



H^ = ^>E, 



WE = EE. 



(6) 



Apart from the bi-orthonormality relation S^^* = I2, it 
holds ^'J'tiyt = H^!^\ EE'^H = fftsst so that one 

rj as metric operator in 



identifies i'^^A 






the pseudo-Hermiticity condition rjH = H^rj [7|]. Ad- 
ditionally to its obvious Hermiticity rj = rf the metric 
can be suitably scaled to be an element of the hyper- 
bolic ("boost") sector of the complex orthogonal group 
50(2,C) i 



I 



r/ 



1 



cos(q:) \y isin(a) 

cosh(/3) ~i sinh(/3) 
i sinh(/3) cosh(/3) 



i sin(a) 
1 



o/3<T„ 



(7) 



with parameter identification sin(a) =: tanh(/3) and 
cos(a) = l/[cosh(/3)]. As final ingredient we fix the no- 
tation for the one-to-one similarity mapping between the 
PT— symmetric Hamiltonian H and its isospectral Her- 
mitian counterpart h = W [Z], El = p~^hp, H^ = php~^, 
p^ = rj, as well as for the unitary eigenvector matrix $ 



/i$ = $£■, 



$t ^$- 



(8) 



The eigenvectors of H and H^ can be regarded as a ^ 
—a, i.e. (3 ^ — /3, mirror symmetrically distorted ver- 
sions of the eigenvectors of the Hamiltonian h 



* = p"i$. 



^ = p$. 



p-\(i)=p{-P). (9) 



The orthogonal initial and final vectors IV'/)! IV'f) in 
the BBJM-brachistochrone model on their turn can be 
considered as eigenstates of a Hermitian spin operator 
Sz — <Jz (a von-Neumann observable with orthogonal 
projector decomposition), whereas the T'T— symmetric 
(non-Hermitian) Hamiltonian H has nonorthogonal 
eigenvectors \E±) and is not a von-Neumann observable. 
Under the one-to-one equivalence mapping [7| from H 
to the Hermitian Hamiltonian h the spin operator Sz 
maps into a non-Hermitian operator s^ — pSzP^^ 7^ s\. 
Hence, the BBJM-brachistochrone system in both repre- 
sentations (H,Sz) and {h,Sz) contains operators which 
are not von-Neumann observables and therefore the sys- 
tem cannot be considered as fundamental. 

Naimark dilation In order to give the BBJM system 
with evolution '0(i) = U{t)ipi a meaning in CQM we 
embed it into a larger purely Hermitian system 



TPit) = U(i)^ 



/; 



^W = 






(10) 



with unitary evolution operator U(i) = [U^'(i)] and 
additional ancilla wave function component x(0- For 
this purpose we construct an auxiliary POVM [5|, [OJ, 
lll| J2k=i^k = h over the nonorthogonal eigenvec- 
tors ([2]), ^ oi H and its adjoint H^ with rank- 
one operators A12 = P\E±{a)){E±{a)\ and ^3^4 = 
f'^\E±{~a)){E±{—a)\. For symmetry reasons all Af^ are 
scaled with the same normalization factor 



/:= 



cos(a) 



1 



y/2 cosh(/3) 



(11) 



Next we note that any rank-one-type POVM built from 



C^ reads in ma- 



a general vector set {xk}k^i : ^k 
trix form Y2ti Xkxl = MM^ = In 
[xi,X2t ■ ■ ^ Xn] € C^^", where M describes a partial ma- 
trix isometry M : C" -^ C^ with n > N. The lifting 
of this isometry to a unitary mapping V : C" —^ C" is 
known as Naimark dilation (extension) [9| and can be ac- 

r M ' 

complished by a suitable embedding M ^^ V = 
L e £^(n-N)xn ^j^j^ additional unitarity constraint 

VVt = VtV = I„. (12) 



For our i 
from eqs. 

where $ 



-<2x4 



setup the matrix M can be read off 



and (| 
$ 
$ 



as M 

= /2< 



$. The natural ansatz 



/ 



P 



P 
X Y 

tion MM^ = p[p'^ + p "i = j^iij + 7] -J = 12 
from dll), dH) and ^) and the constraint ^ fix the 
nonsingular matrices X, Y up to an irrelevant unitary 
rotation a,s X = p, Y ~ —p^^ so that 



$ together with the auxiliary condi- 
+ p"^] = P[ri + rp^] = /2 (following 



V := /K 

= .fWz( 



-P 

-1 



)p 



' p] {h ® *) 



(13) 



The columns of V = [vi, V2,V3, V4] are formed by four or- 
thogonal vectors Vk G 7^4 = C^ which provide the desired 
Naimark dilations of the nonorthogonal eigenvectors of H 
and its adjoint H^ living in 7^2 = C^. Additionally, they 
yield the embedding Ak ^ Pk = Vkvl of the 7^2-POVM 
into the ortho-projector set PjPk = SjkPk, J2k=i ^k — I4 
iiiHA 9\. 

We start the construction of a selfconsistent CQM in 
7^4 by requiring that the original eigenvalue problems for 
H and iJ^ are recovered when the model is restricted 
to the first two rows of V. From relations ([6]) and an 
ansatz /[iJ*, iJ+S] = f[^E,EE] = ME the eigenvalue 
matrix E for the dilated problem can be read off as E :— 

~ ] = I2 ® E. This means that the corresponding 
E J 

dilated Hamiltonian H will have the two eigenvalues E± 

of H and its isospectral adjoint H'' as double degenerate 

eigenvalues. The Hamiltonian H itself can be constructed 

from eq. §i) and HV = VE as H = VEVt so that 



H 



where 



f [I2 <E) {Hf]-^ + rjH) + i(jy ®{H- fft)] 



/2 A -f iuy 



n 



Uq 



A := /^(ffr?-! +ryi/)- ^0/2 + YCOs(a)a, 

This H is Hermitian by construction. In the Hermi- 
tian limit of the original T^T— symmetric Hamiltonian 



H, i.e. for a = /3 = 0, it holds rj = I2 and H re- 
duces to Vi. = I2 ® h — a trivially doubled h. In con- 
trast to the PTQM Hamiltonian H its dilation H remains 
well defined also in the strongly non-Hermitian vanishing- 
passage-time regime ^ where the matrix components of 
H diverge for fixed wq as s — > 00. This regularization 
effect is due to the normalization factor /^ induced in H 
via the auxiliary POVM construction. 

The H— induced unitary evolution in 7^4 is governed by 



the operator U(t) = e 
0, ©, (HID, U{t) = e-**^ 
can be represented as 



itH ^ Yg-^Etyt which 



p ^e **''p and y = uJot/2 



U(i) 



F 
G 



f{l2<»[U{t)rj''+r^U{t)] + 
+i(jy (g) [Uit) - TjU{t)Tj-^] } 

p-^Eot jj^ cos(y) - ia^ sm{y) cos(a)] 
e-'^"* [sm{y) sin(a)cr^] . (14) 



Physically U(i) describes the time evolution of the cou- 
pled brachistochrone-ancilla system (|10p in a Hilbert 
space 7^4 = 7^2 © 7^2 with ■(/'(i) e 7^2 and x{t) G 7i2- 
In order to exactly reproduce the 7^2— evolution ([3|) of 
the BBJM-brachistochrone subsystem 



,-i-Eot 



i^it) = Uit),Pi 



cos{y — a) 
cos{a) y -isin(y) 



(15) 



the initial vector xi € 7^2 of the ancilla subsystem should 
be chosen appropriately. To obtain xi we represent 
^p{t) e 7^4 as 

m = (^ Jfjj ) = e+ ^ ^(i) + e_ ® xit) (16) 



with e 



(1,0)^ and e_ := (0,1)^, define P± 



e± (8) e±. and introduce the projectors P± = P^ (^ I2 
on the brachistochrone (P+) and the ancilla (P-) sub- 
space. The identification rule (fTO|) takes then the form 
'P+'4>(t) = P+U(i)^j = e+ «) -0(0 = e+ (g) C/(t)0/. After 
taking the time derivative one finds from 

P+H^(i) = e+ [A0(t) + ^x{t)] =e+® Ht(;{t) (17) 

a synchronization link[15| between ancilla and brachis- 
tochrone evolution x(i) = n^^{H — A)ip{t) = ?70(i) 
as well as the explicit ancilla evolution x{t) = 

' Initial and final ancilla compo- 



1 



cos(a) \^ _jsin(y _ a) 

nent take then the form xi — vl'i — — j- 



XF = vi^F = j^ 



sin(a) 



W \^isin(a) J' 
with p = e^'^""^ an irrel- 



evant phase factor and ■(/;/ — (1, 0) , tpp = m(Oi 1) 



Discussion For a BBJM-brachistochrone in the van- 
ishing passage time regime ([4]) the ancilla vectors xi 
and xf become colhnear and their common denomina- 
tor cos(a) « £ makes them highly dominant compared 
to ipi,F- This X— dominance remains preserved for the 
normahzed state vector (p := gtp, {4>\4>) ~ 1, with 
g := cos{a)/\/2 and leads to a very small brachistochrone 
component \gt(;{t)\'^ « e^/2 compared to the ancilla com- 
ponent |5x(t)P « 1 — e^/2. As result the geodesic dis- 
tance between the initial and final states in Ti^ becomes 
small (54 = 2arccos(|(^j|^^)|) « 2e. This means that the 
original geodesic distance 62 = 2arccos(|('(/'/|'0F)|) = ti" 
between the initial and final states in the brachistochrone 
subsystem is strongly contracted by embedding the lat- 
ter into the larger Hermitian 7^4— system. Geometrically, 
this follows from the fact that the geodesic distance is 
given by the angle spanned by the corresponding vectors 
on the Bloch sphere [2] and its generalization to higher 
dimensions. In the 7^2— subsystem the vectors t/jj and 
TpF are antipodal and span an angle S2 — tt. Adding a 
much longer vector x orthogonal to ipi and ipp makes 
the resulting cpj w (^/,x)"^ ^md (pp « (V'f,x)"^ al- 
most collinear (^4 ^ in 7^4. In this way the dilated 
model reconciles the Aharonov-Anandan lower bound 
[2, [iSl] on minimal passage times in Hermitian systems 
with the vanishing passage time effect of the BBJM- 
brachistochrone [l| for orthogonal states in the subsys- 
tem. The embedding of the BBJM-system into a higher- 
dimensional Hilbert space can be regarded as a strength- 
ening of the wormhole analogy drawn in [l| for the short- 
ening of the passage time r. A wormhole connection of 
two distant points on a given lower dimensional manifold 
M can be best visualized by embedding M into a higher 
dimensional surrounding A/" D Al so that the correspond- 
ing short distance in Af becomes obvious [U] . 

The representation (fT6|) indicates on a natural interpre- 
tation of the obtained Hermitian system (/)(t) = e~'*^(/)j 
as system of two entangled spin 1/2 particles, i.e. as a 
two-qubit system [5|, with Si = a^^h and S2 = l2®(Tz 
as spin operators of the two spin-subsystems. In order to 
observe the BBJM-brachistochrone effect of the subsys- 
tem one has to prepare an initial entangled state 4>i — 
e+ (g) -0/ + e_ (g) x/j to switch on the interaction Hamil- 
tonian H during the passage time r (assumed as smaller 
than the lower passage time bound t/j = tt/o^o) and to 
evolve 4>j into the final state (pp — e~^'^^4>j. This final 
state has to be analyzed in a two-step measurement. In a 
first (instantaneous) Ei— measurement one selects (filters 
out) the up-component e+ of the first spin. This results in 



a state P+4>p/{P+4>p\T'+4>f 



\i/2 ^ 



e+(^ipp/{ipp\ippy/^ 



and separates the brachistochrone component from the 
ancilla component (connected with the down-component 
e_ of the first spin). In a subsequent E2— rneasurement, 
one analyzes the spin-up and spin-down states of the 
brachistochrone component ■0f/(^f|'0f)^ to recover 



the spin-fiip effect from ^/ = (1,0)-^ to tpp = /iF(0, 1)'^. 

A direct experimental test should be feasible with a 
suitably designed system of entangled photons passing 
an appropriately chosen arrangement of beam splitters, 
phase shifters, and mirrors as implementation of the uni- 
tary operator U(t) = e"*"^^ [lC|- 

Conclusions We have demonstrated that the quan- 
tum brachistochrone for a system with PT— symmetric 
Hamiltonian can be realized as subsystem of a larger 
Hermitian system living in a higher-dimensional Hilbert 
space. The Hermitian system (constructed by Naimark 
dilating an auxiliary positive operator valued measure) 
has the structure of an entangled two-quibit system. 
This opens a way to direct experimental tests on the 
recently hypothesized 'faster than Hermitian' evolution 
in T'T— symmetric quantum systems. 
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